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ABSTRACT 

The ability of new instruments for providing accurate inferences of vector magnetic fields and line-of-sight 
velocities of the solar plasma depends a great deal on the sensitivity to these physical quantities of the spectral 
lines chosen to be measured. Recently, doubts have been raised about visible Stokes profiles to provide a clear 
distinction between weak fields and strong ones filling a small fraction of the observed area. The goal of this 
paper is to give qualitative and quantitative arguments that help in settling the debate since several instruments 
that employ visible lines are either operating or planned for the near future. The sensitivity of the Stokes profiles 
is calculated through the response functions (e.g. Ruiz Cobo & Del Toro Iniesta, 1994). Both theoretical and 
empirical evidences are gathered in favor of the reliability of visible Stokes profiles. The response functions 
are used as well for estimating the uncertainties in the physical quantities due to noise in the observations. A 
useful formula has been derived that takes into account the measurement technique (number of polarization 
measurements, polarimetric efficiencies, number of wavelength samples), the model assumptions (number of 
free parameters, the filling factor), and the radiative transfer (response functions). We conclude that a scenario 
with a weak magnetic field can reasonably be distinguished with visible lines from another with a strong field 
but a similar Stokes V amplitude, provided that the Milne-Eddington approximation is good enough to describe 
the solar atmosphere and that the polarization signal is at least three or four times larger than the typical rms 
noise of 10~^/c reached in the observations. 

Subject headings: Radiative transfer, response functions, magnetic fields, spectropolarimetry, solar magnetism. 



1. INTRODUCTION 

A debate has been maintained for quite long about the ca- 
pability of visible spectral lines to accurately measure weak 
magnetic fields and about whether can they be mistakenly 
confused with strong fields and low filling factors. The dis- 
cussion often reaches levels for the whole community to re- 
think the current means for retrieving the magnetic properties 
of the solar photosphere. Should such confusions be as easy 
and widespread as suggested by several papers, the usage of 
visible spectropolarimeters and magnetographs would advis- 
ably be stopped. As a matter of fact, however, modern visi- 
ble spectropolarime ters, e.g. the Hinod e (Kosu gi et al. 2007; 
iTsuneta et al.l2008b spectropolarimeter dLites et al.l2001l) . are 
providing many new important discoveries whose reliability 
seems not to be easily put in doubt. Here we try to show that 
this is indeed the case and that reasonable accuracy can be 
achieved with visible lines whereby strong and weak fields 
cannot be muddled up easily. 

Advances in infrared spectropolarimetric instrumenta- 
tion have improved significantly our ability to infer the 
magnetic and dynamic properties of the solar atmosphere. 
The large wavelengths and Lande factors of some spectral 
lines make them excellent diagnostic tools because, for 
instance, some may display completely split a compo- 
nents for magnetic fields stronger than, say, 400 G. The 
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advent of infrared data has brought about a controversy 
on the results obtained for internetwork fields when com- 
pared to those coming from visible-line observations. The 
visible Fe I lines at 630 nm have be en indicating a pre- 
dominance of kG fields dSanchez Alm eida & Lites 200^ 
Dommguez Cerdeiia et al.l I2003t [S ocas-Navarro & LitesI 
2004), whi l st the infrared Unes at 1565 n m suggested hG 
fields (Hn T995"; 'Lin & Rimmele' 'Collados ' '2001 

Khomenko et al. 2003; Dommguez Cerden a et al. 200™. 
Whereas a decade ago IWestendorp Plaza et al.l ( il998h 
provided evidence on the accuracy of Milne-Eddington 
(ME) inversions of the pair of Fe I lines at 630 nm for 
weak fields, some recent works even raise serious doubts 
about the reliability of visible-line inferences through in- 
ver sion dMartmez Gonzalez et al. 200 0). T he suggestion 
by ISocas-Navarro & Sanchez Almeida ( i2003l) about a bias 
of visible Unes to strong fields and infrared lines to weak 
fields is too disappointing because it somehow admits an 
inability of visible lines to diagnose an important fraction 
of the solar magnetism. There certainly are differences in 
the sensitivities of lines in both wavelength regions but this 
should not mean that we must abandon spectropolarimetric 
measurements in the visible. There is much more information 
available if simultaneous observations in the visible a nd the 
infrar ed are used (e.g., i Cabrera Solana et al. 2 0061 l2007l 
I2OO8I; iMartmez Gonzalez et al.. ,2008; Beck ,200S|k But is 
there any hope that visible Une observations can provide 
reliable inferences when weak fields are observed? A solu- 
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tion to this controversial dispute is urgently needed because 
infrared wavelengths are not available to many polarimeters. 
Indeed some new instruments have been devised to use 
visible lines, namely, the Hehose ismic Magnetic Imager 
(IScherrer & SDO/HMI TeamI 12002 ) aboard Solar Dynam- 
ics Observatory, the Imag ing Magnetograp h experiment 
dMartmez Pillet et al] |2004 aboard Sunrise dBarthol et alj 
120081) . and the Polarimetric and Helioseismic Imager 
aboard Solar Orbiter. The results by lOrozco Suarez et aLl 
(|2007a'b"c) from a thorough analysis of magnetohydrody- 
namic simulations and real data have shed some light to the 
problem because their results indeed reconcile visible and 
infrared measurements. They concluded on the reliability of 
visible-line results if enough spatial resolution is available 
and stray light from the surroundings is appropriately taken 
into account. 

The present paper tries to provide additional arguments to 
help discriminate whether or not visible spectropolarimetry is 
able to disentangle the two extreme cases when inferences are 
made through inversion techniques of the radiative transfer 
equation. In Sect. |2]we clarify a number of points regarding 
inversion techniques and the underlying model atmosphere 
and physical assumptions; in Sect. [3] we discuss on the in- 
trinsic ability of visible lines to distinguish weak from strong 
fields; in Sect. |4]we calculate the uncertainties of regular ME 
inversions if only noise is a limiting factor. These uncertain- 
ties can be considered as detection thresholds for the corre- 
sponding model parameters. Finally, in Sect.|5]we summarize 
the conclusions of this work. 

2. MODELS, SCENARIOS, AND INVERSIONS 

The vocabulary related to inversion techniques, and to as- 
trophysical measurements in general, is sometimes mislead- 
ing to some extent. We often speak about model atmospheres, 
about scenarios, about hypotheses without clarifying their 
meanings. Let us, therefore, specify what we mean in this 
paper for each term and how we use all the concepts involved. 
On the other hand, practice teaches that failures occur from 
time to time and inversions might eventually not be satisfac- 
tory. Why does an inversion code occasionally fail to find a 
reliable solution? Is there anything intrinsic to the procedure 
that disable the success capabilities? We address these two 
issues in this section 

First, a distinction should be made between a model atmo- 
sphere, that is, the set of parameters that fully characterize 
the solar medium from which the Stokes spectrum is coming, 
and the physical scenario, i.e., the set of physical assumptions 
about the composition of such model atmosphere and hence 
the means for defining which model parameters are relevant. 
For example, a physical scenario can be such that the physical 
quantities are assumed to be constant with depth; if we in- 
clude the effects of stray light or of a partial filling of the res- 
olution element, the corresponding model atmosphere can be 
the set of nine ME parameters' plus the filling fraction of the 
magnetic component. A given set of Stokes profiles could be 
inverted after assuming two different scenarios; for instance, 
a second scenario would be that the magnetic and velocity 
field vary with depth through the atmosphere. Should both 
of them produce equally good fits to the observations (similar 

' So and Si for the constant term and the slope of the source function, rio 
for the Hne-to-continuum absolution coefficient ratio, AAd for the Doppler 
width of the line, a for the damping parameter of the line, B, 7, and i/p for the 
strength, inclination, and azimuth of the vector magnetic field, v^os for the 
line-of-sight velocity, and a for the magnetic filling factor 



values of the merit function), the Occam's razor principle is 
usually applied: the simpler, the better. In other words, prac- 
tice recommends to be careful with involved scenarios (e.g. 
gradients along the line of sight or several components within 
the resolution element) whenever a simple one (e.g. ME) pro- 
duces good fits to the observations. Of course, fit quality is 
a subjective concept but can be quantified in terms of (x^) 
merit-function values: we can set a threshold for the differ- 
ences between observed and synthetic (best-fit) profiles; if 
reaches smaller values than this threshold, we can say we are 
satisfied with the inversion. Let us explain our point other- 
wise: imagine that the second scenario does not lead to a clear 
minimum of the merit function because two different sets of 
model parameters yield no significantly different x^- In such 
a case, we cannot say that the observables (e.g., the visible 
line Stokes profiles) do not provide reliable diagnostics unless 
we try with simpler scenarios like the first one with constant 
physical quantities. It may well be that this scenario provides 
good fits to the data.^ In Sect. I3.2l we show one such case. 

Second, there are only two reasons for a Levenberg- 
Marquardt inversion code to fail, namely, that the merit func- 
tion has no well-defined global minimum in the hyperspace of 
parameters or that the updating strategy for the model atmo- 
sphere-' is not well suited to the problem. Let us discuss both 
reasons separately. 

Consider two or more model atmospheres within the same 
scenario. If they produce equally good fits (i.e., equally low 
values of the merit function), then either the observational 
noise is such that hides the true minimum, or the underly- 
ing physical scenario is not appropriate to describe the so- 
lar atmosphere, or both. We take for granted that noise ham- 
pers the measurements. Indeed it governs the threshold above 
which the model parameters can be reliably inferred (see Sect, 
m. Discarding a situation where the noise is too high, how- 
ever, let us speak about the scenario. By ambiguity of the 
scenario we mean that the physical constraints behind it can 
be such that the number of free parameters is too high for 
the information content available from the observables or that 
some parameters cannot be reli ably retrieved. For instance, 
iMartmez Gonzalez et al. I (l2006 l) have shown that the pair of 
Fe I lines at 630 nm is not able to provide a single model for 
a scenario in which two depth-dependent atmospheres, one 
magnetic and another non magnetic, fill a spatial resolution 
element of about 1". This can be an example of the first case: 
several combinations of the atmospheric quantities yield the 
same value. As for the second case, a good example is 
found in the well-kno wn trade-offs among 770, AAj , and a 
of ME inversions (see IWestendorp Plaza et al1ll998l Orozco 
Suarez & Del Toro Iniesta 2007): several sets of such three 
parameters may give fits with the same quality without chang- 
ing the magnetic and velocity parameters. This means that the 
thermodynamic parameters of a ME atmosphere are not reli- 
able while the magnetic and dynamic are. 

Let us suppose now that the physical scenario explains the 
observations, that is, that a clear global minimum should 
exist. To reach it, the usual strategy is to increase or de- 

^ By good fits we understand both values below the established thresh- 
old and solution uniqueness or independence of the initial guess. 

' We hereafter restrict the discussion to Levenberg-Marquardt algorithms 
as are among the most used ones in practice. The fudge parameter A is the 
technical means for these algorithms to update the model atmosphere by ap- 
propriately combining the advantages of the inverse-Hessian and the steepest- 
descent methods. The latter is in general better when far from the global 
minimum and the second when close to it. 
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crease the A parameter by a factor 10, depending on the 
divergence or convergence of the algorithm. If the code 
reaches a local minimum, then the fudge parameter has to 
be large enough to favor bigger perturbation steps and over- 
come that minimum. Sometimes, however, the updating strat- 
egy is such that the fudge parameter is unable to get rid 
of the local minimum. For example, in the MILOS pro- 
gram dOrozco Suarez & Del Toro Iniesta Il2d07h A monoton- 
ically decrease until it starts oscillating and gets trapped. Al- 
ternative updating strategies as well as initial values for the 
fudge parameter can be devised to increase the efficiency of 
convergence, but ultimately they are somehow based on the 
user's experience. The only automatic option is to stop the 
inversion procedure as soon as oscillations appear. 

No major technical problem seems to exist, then, provided 
that an appropriate scenario and a good fudge-parameter strat- 
egy have been selected for inversion codes to reach good fits 
to the observations. The distinction between different model 
atmospheres within a single physical scenario (e.g. the dis- 
tinction between weak and strong magnetic fields) only de- 
pends on the information content in the Stokes profiles. Such 
a content depends in its turn on the spatial resolution (filling 
fraction of pixel occupation) and on the noise of the observa- 
tions. The conclusions, hence, cannot be general for given 
spectral lines but specific for a given observation. For in- 
stance, Bellot Rubio & Collados (2003) distinguish between 
weak and strong fields provided the observational differences 
are significantly above the noise. 

3. DETECTABLE SIGNATURES OF WEAK AND STRONG FIELDS 

This section is devoted to show the differences between 
Stokes profiles formed in atmospheres with weak magnetic 
fields plus large filling factors and strong magnetic fields but 
small filling factors. Such differences enable the inversion 
codes to reliably distinguish between the two cases as we are 
going to see through a double approach. First, we use syn- 
thetic profiles and theoretical arguments. Second, a sample 
case with actually observed Stokes profiles is shown. They 
have been inverted with the MILOS code. This code deals 
with individual lines but, like other ME inversion codes ex- 
isting in the literature (e.g. Skumanich & Lites 1987), MI- 
LOS can invert simultaneously specific pairs of lines. This 
is the case of the two Fe I lines at 630 nm. On the one 
hand, inversion of lines not being normal Zeeman triplets like 
that at 630.15 nm is possible by computing the full Zeeman 
pattern and considering the components altogether, following 
lAuer et all (11977 ). On the other hand, the two lines belong to 
the same multiplet and, hence, are fairly similar in everything 
except for their line strength. This property allows the user 
not to increase the number of free thermodynamic parameters 
by simply consideri ng the ratio of their oscillator strengths 
(ISkumanich & Litesl[l987h and, thus, of their ryo's: just one 
rjo, one a, one AAd, one Sq, and one are retrieved con- 
sistently with one B, one 7, one (p, and one a. In any other 
technical aspect, MILOS is a regular Levenberg-Marquardt 
non-linear, least-squares algorithm that inverts the full radia- 
tive transfer equation under the assumption of the ME approx- 
imation using measurements of the four Stokes parameters to 
better constrain the inversion."* 

3.1. The synthetic case 

It is programmed in IDL language and is publicly available on our Web- 
site, http://spg.iaa.es/download.asp 



Consider two ME atmospheres whose thermodynamic pa- 
rameters are the same. Let us take for instance those coming 
from a fit to the intensity profile of the Fe I line at 630.25 
nm i n the Fourier Transform Spectrometer A tlas of the quiet 
Sun (Brault & NeckelllI98"7t lNeckell ll999Y Two different 
atmospheres will be built by assuming a different magnetic 
field strength: 200 and 1500 G. The magnetic inclination and 
azimuth, 7 and <j), are taken equal to 45° for both model at- 
mospheres. The specific values have nothing to do with the 
qualitative results and, in fact, intend to be representative of 
the general problem. From these two atmospheres different 
Stokes profiles emerge as can be seen in Fig. [T] Solid and 
dotted lines correspond to the weak and strong field cases, re- 
spectively. The shapes of the profiles are clearly dissimilar 
and the Zeeman splitting is more conspicuous for the strong- 
field profiles. A filling factor unity is assumed in both cases. 
The differences are so significant that distinguishing between 
the two cases is fairly easy. Inversion techniques should have 
no problem in inferring either of the two models. 
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Fig. 1. — Emergent Stokes profiles from three different ME atmospheres 
with the same thermodynamic parameters. Differences lie on their magnetic 
field strengths and filling factors. Check line coding in the figure inset (right, 
bottom panel). 

It could be argued that a filling factor a = 0.133 in the 
strong-field case would be equivalent to the weak-field at- 
mosphere because the magnetic flux is the same in the two 
models. However, this would disregard the enormous abso- 
lute differences in the profile shapes: as can also be seen in 
Fig.[Tl they can be as large as 1% in Stokes Q/Ic, as 2% in 
Stokes U/Ic, and as 8-10% in Stokes I/I^ and V/Ic- Hence, 
the differences are between one to two orders of magnitude 
larger than the typical noise (I0~^lc) achievable by modern 
spectropolarimeters and are thus discernible by current in- 
version codes. Of course, our conclusion depends on the 
validity of the scenario as the belief on a better accuracy 
for the flux determination than for the field strength deter- 
mination roots in the proportionality between magnetic flux 
and Stokes V. This proportionality is true in a very narrow 
range of field strengths, however: we know, indeed, that strict 
proportionality only takes place while second-order effects 
are negligible; but the mere existence of linear polarization 
breaks this condition: Stokes Q and U are of second order 
jLandi degrinnocenti1 ll992h . Both circular and linear po- 
larization scale with a in the weak field regime but while V 
scales with B, Q and U go as B^ and are even more sensitive to 
changes in the magnetic field strength. This means that wher- 



4 



Del Toro Iniesta, Orozco Suarez, & Bellot Rubio 



S = 200 G 



-200 -100 100 200 

X (mi) 





l.^(10m/s)-' 




B (10 G)"' 




y (■)-' 




200 


-100 100 200 
A (mA) 



10 



-10 : - 



-200 -100 100 200 -200 -100 100 200 

X (mi) A (mJ) 




-mh ^ ^ ^ d -aIl . . ^ c 

-200 -100 100 200 -200 -100 100 200 

\ (mi) i (mi) 



Fig. 2. — Response functions of the Fe I Stokes profiles at 630.25 nm 
(normalized to the continuum intensity) for the weak- (four top panels) and 
strong-field (four bottom panels) ME models versus wavelength. RFs to vl^s 
perturbations (solid, black lines), to B perturbations (blue, dotted lines), and 
to 7 (red, dashed lines) are presented. Perturbations of 10 ms"', of 10 G, and 
of 1 ° are assumed. 

ever we measure linear polarization, the proportionality is not 
exact. Therefore, second-order effects are enough for getting 
accurate inferences provided they are bigger than the noise. 
This is the case, for example, of the weak fields detected in 
the inter network by the spectropolarimeter aboard the Hinod e 
satelUte dLites et al.ll2007ll2008HOrozco Suarez etani2007bl) . 
that had hitherto gone unnoticed from the ground. Even in 
the case of strictly longitudinal fields (i.e. with no linear 
polarization), proportionality breaks down soon. Stokes V 
ceases to scale with field strength at fairly moderate strengths: 
a linear and a parabolic fit to the Stokes V amplitude as a 
function of the (longitudinal) field strength differ by 1-2% of 
the continuum intensity between and 400 G. For 500 G, 
the differences reach a 6%.^ All these differences are sig- 
nificantly larger than the typical noise of 10"^ 7^. Hence, as 
long as we have several wavelength samples (i.e. unless we 
perform single-point magnetometry), the chances of distin- 
guishing non linearities are quite high. Numerical tests of the 
ability of ME inversions of the Fe I lines at 630 nm have al- 
ready been carried out by Westendorp Plaza et al.l (^998) with 
synthetic profile s in both constant and depth-d ependent atmo- 
spheres and by lOrozco Suarez et aLl (l2007ah with synthetic 
profiles from numerical MHD simulations of the photosphere. 
We further check and confirm the reliability of such inver- 
sions with the Monte-Carlo simulations we report in the Ap- 
pendix. There, a total of 10^ ME profiles have been inverted 

^ The quoted differences correspond to the following ME parameters fixed: 
Tja = 22, AAo = 30 mA, a = 0.03, and v^os = 7 = = 0. The values are 
typical for visible Unes like those of Fe I at 630 nm. Small changes in these 
parameters give very similar results. 



after adding them random noise with a rms value of 10"^ Ic- 

The most important feature from the measurement point of 
view, however, is t he sensitivities of the profi les as given by 
response functions dRuiz Cobo & Del Toro In iesta 1994). Re- 
sponse functions (RFs) of the Fe I line at 630.25 nm to per- 
turbations of the LoS velocity (solid lines), of the magnetic 
field strength (dotted lines), and of the magnetic inclination 
(dashed lines) are plotted versus wavelength in Fig. |2] The 
RFs for the weak-field model are in the four top panels and for 
the strong-field model in the four bottom panels. Note that the 
scales are the same for each Stokes parameter in the two mod- 
els. Units have been selected to represent perturbations of 10 
ms"', of 10 G, and of 1°, respectively, since the response to 
a 1 "perturbation can be as large as (or even larger than) those 
to 1 G and lms~'. The spectral line responds differently in 
the two models to the same perturbation in the atmospheric 
parameters. Moreover, the various wavelengths react differ- 
ently to the same perturbation in the two atmospheres: both 
the shape and the amplitudes of the RFs are different. This 
is the key that allows inversion codes to discriminate between 
the two models. 

The relative accuracies of the inferences in the two cases 
are determined by their relative sensitivity or RF values. The 
difference in maximum sensitivities, however, is small: as 
seen in Fig. |2] 3 G can be as significant in one of the mod- 
els as 1 G is in the other. The RF amplitude is even larger for 
Stokes Q and V in the weak model than in the strong model. 
Moreover, it is not only the maximum sensitivity what mat- 
ters but the effect of the whole profile. Therefore, although 
the common belief is that strong fields are better retrieved 
than weak ones, inferences in the weak-field case can even- 
tually be more accurate if the fields are not very horizontal 
as we are going to see in the next section. The point we 
want to stress in here is that the neat shape differences of the 
RFs help inversion codes disentangle strong from weak mag- 
netic fields.^ This is a similar situation to that giving rise to 
a discrimination between the magnetic and dynamic param- 
eters from the thermodynamic parameter s in ME inversions 
(lOrozco Suarez & Del Toro Iniesta Il2()07l) . 
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Fig. 3. — Stokes / (left) and V (right) RFs to v'loS (solid, black fines), to 
B (dotted, blue lines), and to a (dashed, red lines) in the weak-field case. 
Perturbations of 10 ms^'for vl^s. of 10 for B, and of 0.1 for a have been 
assumed. 

The fact that filling factor and magnetic field strength influ- 
ence differently the Stokes profiles can also be assessed with 
RFs. As shown in Fig. [3] in the weak field model the sen- 
sitivity of Stokes y to a perturbations has exactly the same 
shape than the sensitivity to B perturbations (the same occurs 
for Stokes Q and U). This is so because, at these strengths. 

Inversion codes like SIR <Ruiz Cobo & Del Toro Iniestalll993) or MI- 
LOS are based on response functions. Other codes not explicitly programmed 
making use of RFs do indeed use them imphcitly (e.g.. lDel Toro Iniestal2003l) 
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V scales with both a and B. (Both RFs have a V-shaped pro- 
file.) This is not the case, however, for Stokes /: the RFs of 
Stokes / to a and B are very different. In fact, while the re- 
sponse to a perturbations keeps a V-like shape, the RF to B 
perturbations has three lobes. ^ Even in the weak field case, 
then, thanks to the Stokes / behavior, the effects of a and B 
can be discerned. Therefore, the reliability of the separation 
be tween the two cases r oots in Stokes /, as first pointed out 
bv lOrozco Suarez et al.l (l2007d) . The intensity cannot be fit- 
ted by modifying one of the two parameters regardless of the 
other 
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Fig. 4. — Stokes / and V profiles of the Fe I lines at 630 nm as observed 
with the spectropolarimeter aboard Hinode in a network (upper panels) and 
an internetwork (middle panels) points. Observations (solid lines) and fits 
(dashed lines) are displayed. Their differences are also shown at the bottom of 
each panel. The bottom panels show the difference (network - internetwork) 
between the profiles of the two points. Note that the internetwork polarity has 
been artificially changed for calculating this difference. 



' Since the observed spectrum is a linear combination of the magnetic and 
non-magnetic atmospheres, the RF to a perturbations is simply the differ- 
ence between the magnetic and non-mag netic Stokes profiles, as shown by 
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3.2. Observed Stokes profiles: the real case 

Problems of distinction between the two cases invariably 
occur in real observations when polarization signals are weak. 
As an illustration, we have chosen the Fe I lines at 630 nm 
as observed in two points, one from the network and another 
from the internetwork. The data come from the Hinode oh- 
servations already u s ed by [Lites et all |2007, 2008) and by 
lOrozco Suarez et aH (l2007 c b). The circular polarization is 
less than or equal to 1% of the continuum intensity in both 
points. The linear polarization is buried within the noise. The 
Stokes profiles of both points are shown in Fig. |4] (four first 
rows) along with their differences (bottom panels). The ob- 
served (solid line) profiles and the fits (dashed line) are in- 
cluded. The quality of the fits is very good as indicated by 
the residuals: those for the V profiles hardly exceed the 3cr 
level marked by the horizontal dotted lines, (cr = 10~^/c is 
the typical noise level of the Hinode spectropolarimetric ob- 
servations.) The ME inversions carried out with the MILOS 
program give (160 ±2)° and (50 ±2)° for the field inclina- 
tion, 7; 0.1 30 ±0.005 and 0.25 ±0.02 for the magnetic filling 
factor, a; (2.8 ±0.1) kms-'and (2.0±0.1) kms-'for the line- 
of-sight velocity, vlos; and (1400 ±200) G and (280 ± 80) G 
for the magnetic field strength, B. The uncertainties are ob- 
tained after a Montecarlo experiment repeating the inversion 
1000 times and are indeed consistent with the analytic estima- 
tion of uncertainties that we describe in the following section. 

Although the amplitudes are fairly similar (their differences 
are larger than 3(t in just a few wavelength samples), regard- 
less of the magnetic polarity, a net difference in the width of 
the Stokes V lobes can easily be seen. Therefore, the two 
Stokes V profiles are by no means proportional and the lin- 
earity hypothesis with the magnetic flux is not valid at all. 
In fact, the network profile corresponds to a flux four times 
larger than that of the internetwork profile. Nevertheless, al- 
though the subtle lack of proportionality should be detectable 
by the code, the conspicuous difference in Stokes / (well 
above the 3f7 level) points to very different field strengths. 
Since the field inclination angles are almost supplementary, 
the observed / difference can only come from differences in 
vloS, oi, B, and/or the thermodynamics. According to the in- 
version, the network point is shifted to the red with respect to 
the internetwork one. Hence, the expected asymmetric differ- 
ence should be of opposite sign to that observed. The thermo- 
dynamic parameters are known to cross-talk among them (e.g. 
iOrozco Suarez & Del Toro Iniesta 1 120071) and, hence, their 
differences cannot give account of the profile differences. 
Therefore, most of them should come from B and a. In the 
previous section we have provided analytic evidence about 
neatly different sensitivities to B and a perturbations but let 
us for the moment assume that this is not the case and that 
B and a could have been fully confused: if a had to be 0.25 
for both profiles, the network field strength would turn out to 
be 700 G, still very different to the internetwork value, and 
certainly not in the weak field regime. Finally, a further ev- 
idence is the Stokes V amplitude ratio between the two lines 
of the Fe I multiplet. Such a ratio suggests a strong field for 
the network profile and a weak field for the internetwork one. 

4. UNCERTAINTIES 

As is well known, and already mentioned in Sect. |2l noise 
can be a limiting factor to the accuracy of the results of any 
measurement technique and, in particular, of inversion codes. 
In this section we provide quantitative estimates of the un- 
certainties in the retrieved physical quantities induced by the 
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noise if the observations are analyzed through a ME inversion 
code that assumes a magnetic atmosphere filling a fraction a 
of the resolution element. In such a case, the observed Stokes 
spectrum is classically expressed as 

lobs = aim + (l-a)Inm, (1) 

where Im is the Stokes vector emerging from the magnetic 
atmosphere and Inm is obtained as a given average of the pixel 

surroundings^ 

Following ISanchez Almeidal (Il997h . assume that all uncer- 
tainties in the m physical parameters contribute in a similar 
amount to the final noise. (Actually, our assumption is that 
noise imparts equal, Gaussian-distributed uncertainties to the 
m parameters that are sought.) In such a case, the variance of 
the j wavelength sample in the i-th Stokes parameter** can be 
written as^ 

alj=maHRljf<Tl (2) 

where a stands for the magnetic filling factor, RJ j is the RF 
of Stokes / at wavelength j to perturbations of the model pa- 
rameter X, and crj is the variance of that parameter. Summing 
up for all Stokes parameters and wavelengths, Eq.|2]becomes 

tt-f.=-^-'-^tt(Kjf, (3) 

i=i j=i i=i j=i 

where «a is the number of wavelength samples. 

According to Del Toro Iniesta & CoUadosI (l200Cll) . if all the 
modulated measurements (rip) needed to derive the Stokes pa- 
rameters have the same variance, a^, due, for instance, to pho- 
ton noise, then 

o',w= — ^,Vj = l,...,nA, (4) 

rip ef 

where e,- is the polarimetric efficiency for the i-th Stokes pa- 
rameter. 

Using Eq. (|4|i, the uncertainty of the x parameter can finally 

be cast as 

o-.( = , = ■ (5) 

The above formula gives an estimate for the noise-induced, 
i.e. random, effects and no systematic errors are included. 
The latter can b e larger as already show n through numerical 
experiments by iWestendorp Plaza et al.l ( 1998). Our formula, 
however, illustrates very well how the noise on the (modu- 
lated) polarization measurement influences directly the accu- 
racy in any inferred parameter. Obviously, the better the po- 
larimetric efficiencies of the instrument, the smaller the uncer- 
tainties. That is also the case for RFs: the larger the RFs (the 
sensitivity) the smaller a^. Values from Eq. (|5]l are plotted in 
Fig. |5]for a range of magnetic field strengths (from through 
1850 G) and magnetic inclinations (15° in dotted lines, 45° in 
solid lines, and 75° in dashed lines). The model atmosphere is 
an average of those obtained aft er the inversion of more tha n 
two million quiet-Sun profiles (lOrozco Suarez et aDl2007bh . 

Index (' runs from 1 through 4, corresponding to Stokes /, Q, U, and V, 
respectively. 

' Within the ME approximation, RFs are direct partial derivatives 
of the Stokes spectrum vyith respect to the model parameters (see 
lOrozco Suarez & Del Toro Iniesta 112007ft . Obviously, only Im depends on 
the free parameters. 
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Fig. 5. — Uncertainties for the retrieved field strength (top left panel), 
for the field inclination (top right panel), for the LoS velocity (bottom left 
panel), and for the filling factor (bottom right panel), according to Eq. js). 
as functions of the magnetic field strength. Line colors and types indicate 
different field inclinations; see the inset. 

nx = 9l, Up = 4, ej = 1, 1/3, 1/3, 1/3, a = 1, and a = 10" V^. 
The best accuracies are obtained for vertical fields in all the 
strength regimes, except for the inclination that is (naturally) 
better determined when the field is close to the horizontal. 
The improvement on the velocity determinations for fields 
stronger than approximately 900 G for t he 7 = 45°, 75° cases 
is due to a well-known shape effect (e.g. lCabrera Solana et al.l 
2005): with these inclined magnetic fields, the tt component 
of the line starts to become prominent on Stokes /, making it 
deeper and narrower That is also the reason for the overall 
increase in the vloS uncertainties with the field strength: the 
lines become broader as the strength grows. The neat decrease 
of aB with B in the 7 = 45°, 75° cases is due to the appearance 
of Stokes Q and U signals. The slight but appreciable in- 
crease of the uncertainty with B in the 7=15° case might be 
due to a net decrease in information: when the field is strong, 
larger strengths mean larger separation of the V lobes but the 
variation in amplitude is small. Hence, the field strength un- 
certainties in the three inclination cases converge for strong 
fields. The behavior of the inclination and filling factors is 
supposed to be more natural: determinations are better when 
B is large. Nevertheless, as we advanced in the previous sec- 
tion, it is the effect of all the four Stokes profiles that is rele- 
vant to the final inference; arguments based on just one Stokes 
parameter may fail. It is finally noteworthy that the accuracy 
in each parameter is inversely proportional to the magnetic 
filling factor, according to Eq. (|5]l. Hence, the ordinate scale 
of the figure should be multiplied by 10 if a = 0.1. In such 
a case, and even if our estimates were wrong by a 100%, the 
expected uncertainties support our conclusion that strong and 
weak fields can be distinguished with visible lines provided 
that a ME atmosphere partially filling the resolution element 
can be assumed as a model of the solar photosphere that is 
observed with a noise of the order of 10~''/c. 

5. CONCLUSIONS 

In this paper we have gathered evidence in favor of the abil- 
ity of visible Stokes profiles to provide reliable information 
about the magnetic and dynamic properties of the solar pho- 
tosphere and, in particular, to distinguish between strong and 
weak magnetic fields with similar Stokes V amplitudes. Our 
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results complement those bv lOrozco Suarez et all (l2007allblB 
and agree w ith them. Most of the problem s found by other 
authors (e.g. lMartmez Gonzalez et aini2006h can be ascribed 
directly to using scenarios that are too complex for the in- 
terpretation of the observations. The additional information 
provided by infrared lines or by the combined use of visible 
and infrared lines enables us to attack the problem by assum- 
ing complicated scenarios with several components within the 
resolution element whose physical quantities vary with depth 
and so on. These scenarios may not be appropriate for visible 
lines or for observations with limited wavelength sampling. 
Such observations should be analyzed in terms of less ambi- 
tious (hence simpler) scenarios where, probably, the number 
of atmospheric components is reduced to one plus a stray- 
light contribution and where physical quantities are constant 
with depth. 

We have shown that linearity between Stokes V and a and 
B cannot be any longer an excuse for people not trusting the 
results based on the inversion of visible lines. The reason is 
four-fold: linearity between V and aBcos7 takes place in a 
very narrow range of values and cannot be taken for granted; 
circular polarization scales with B whilst linear polarization 
does it with B-\ linear polarization helps a lot for measuring 
B\ the sensitivities of Stokes I to B and a are conspicuously 
different in any field regime and help to disentangle the ef- 
fects of both parameters; and the uncertainties in the model 
parameters induced by the noise of modern spectropolarime- 
ters set detection thresholds that are small enough. (These 
theoretical thresholds are fully compatible with the numerical 
results from Monte-Carlo simulations.) Occasional mistakes 
(e.g., disproportionally large field strengths with extremely 
small filling factors) obtained with existing inversion codes 
cannot then be ascribed to intrinsic difficulties based on ra- 
diative transfer theory but likely to either specific numerical 
and/or programming problems, or to ambitious analyses with 
too many free parameters. 

APPENDIX 

To further check the reliability of our ME inversions of the 
Fe I pair of lines at 630 nm, in this appendix we perform the 
following numerical experiment: first, we generate a refer- 
ence set of Stokes profiles; second, we add Gaussian noise 
with a rms value of IQ~^ 1^; and, finally, we invert the Stokes 
profiles with the MILOS code. According to Eq. [T] a mag- 
netic atmosphere occupies a fraction a of the resolution ele- 
ment and a non-magnetized one occupies the rest. As for the 
magnetic component, 100000 ME model atmospheres with a 
uniform random distribution of vector magnetic fields (field 
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Strengths from to 3000 G and inclinations and azimuths 
from to 180°) have bee used. The remaining model pa- 
rameters are those coming from the inversion of the FTS at- 
las data as in the main text (SecL lTTT i. The filling factor has 
been considered with a uniform random distribution as well, 
with values from a = to 1 . The non-magnetized atmosphere 
is assumed to be known a priory. Same thermodynamic pa- 
rameters as those for the magnetic component have been as- 
sumed. An extra broadening to the profile using a macrotur- 
bulent velocity of 1 km s~' has been used which is consistent 
with a contribution mostly coming from the surroundings of 
each pixel. 

All inversions have been carried out with the following ini- 
tialization: Ao = 1 (fudge parameter), = 0.02, 5i = 1, 770 = 4, 
B = 800 G, 7 = 45°, X = 45°, AAd = 26 mA, vlos = km s"', 
a = 0.15, and a = 0.6. We have used a maximum of 30 itera- 
tions per pixel. As is usual in practice, only profiles displaying 
at least one of the polarization signals larger than 3 10"^ 4 are 
inverted: 5% of the profiles have hence been discarded. 

The upper panels of Fig. |6] show the difference between the 
inferred parameters and the reference model. Each point cor- 
responds to an individual trial of each of the 100000 ME in- 
versions. Over-plotted are the corresponding mean (solid line) 
and rms values (dashed line). The panels show the values for 
the magnetic field strength, its inclination, azimuth, and the 
filling factor The results indicate that the inferences are ac- 
curate enough for all the four model parameters. For fields 
of 100 and 200 G, the mean errors are 30 G and -3 G for the 
field strength, 1 % and 0.9% for the filling factor, 3° and 2° 
for the field inclination, and -40° and -3° for the field az- 
imuth, respectively. The corresponding rms values are: 30 G 
and 10 G, 2.2 % and 1.7%, 25° and 12°, and 45° and 40°. The 
larger values for the azimuth errors are the result of the very 
low linear polarization signal produced by these very weak 
fields. Moreover, the shapes of the various distributions are 
fully compatible to those in Fig. |5] calculated from Eq. (|5]l. 
Obviously, since we have here used a broader range of values 
for the various model parameters, the scatter of the different 
panels tends to be slightly larger than that for the theoretical 
predictions. 
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